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1 Introduction 

Calculus of variations on time scales (we refer the reader to Section 2 for a 
brief introduction to time scales) has been introduced in 2004 in the papers by 
Bohner [2] and Hilscher and Zeidan [4] , and seems to have many opportunities 
for application in economics [1] . In both works of Bohner and Hilscher&Zeidan, 
the Euler-Lagrange equation for the fundamental problem of the calculus of 
variations on time scales, 

,b 

C[y(-)]= L(t,y°(t),y A (t)) At ^ mm, y(a) = y a , y(b) = y b , (1) 

J a 

is obtained (in [4] for a bigger class of admissible functions and for problems 
with more general endpoint conditions). Here we generalize the previously 
obtained Euler-Lagrange equation for variational problems involving delta 
derivatives of more than the first order, i.e. for higher-order problems. 
We consider the following extension to problem (1): 
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C[y(.)} = / L(t, y° (t), y° " A (t), ...,y° A ~ (t), y A (t))At min, 

J a 

y(a) =y a , y(p r_1 (6)) =y b , 
y Ar -\a)=y:-\ y A ^ (^(6)) = yl~\ 

(P) 

where y a A (t) G M™, i G {0, . . . , r}, n, r G N, and £ belongs to a time scale 
T. Assumptions on the time scale T are stated in Section 2; the conditions 
imposed on the Lagrangian L and on the admissible functions y are specified 
in Section 3. For r = 1 problem (P) is reduced to (1); for T = K wc get the 
classical problem of the calculus of variations with higher-order derivatives. 

While in the classical context of the calculus of variations, i.e. when T = R, 
it is trivial to obtain the Euler-Lagrange necessary optimality condition for 
problem (P) as soon as we know how to do it for (1), this is not the case on 
the time scale setting. The Euler-Lagrange equation obtained in [2, 4] for (1) 
follow the classical proof, substituting the usual integration by parts formula 
by integration by parts for the delta integral (Lemma 2). Here wc general- 
ize the proof of [2, 4] to the higher-order case by successively applying the 
delta-integration by parts and thus obtaining a more general delta-differential 
Euler-Lagrange equation. It is worth to mention that such a generalization 
poses serious technical difficulties and that the obtained necessary optimality 
condition is not true on a general time scale, being necessary some restric- 
tions on T. Proving an Euler-Lagrange necessary optimality condition for a 
completely arbitrary time scale T is a deep and difficult open question. 

The paper is organized as follows: in Section 2 a brief introduction to 
the calculus of time scales is given and some assumptions and basic results 
provided. Then, under the assumed hypotheses on the time scale T, we ob- 
tain in Section 3 the intended higher-order delta-differential Euler-Lagrange 
equation. 

2 Basic Definitions and Results on Time Scales 

A nonempty closed subset of R is called a time scale and it is denoted by T. 
The forward jump operator a : T — > T is defined by 

a(t) = inf {s G T : s > t}, for all t G T, 

while the backward jump operator p : T — > T is defined by 

p(t) = sup {s G T : s < t}, for all t G T, 

with inf = sup T (i.e. a(M) — M if T has a maximum M) and sup = inf T 
(i.e. p(m) = m if T has a minimum to). 
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A point t G T is called right-dense, right- scattered, left-dense and left- 
scattered if <j(t) = t, er(t) > t, p(t) — t and p{t) < t, respectively. 

Throughout the paper we let T = [a, b] n TT with a < b and T a time 
scale containing a and b. 

Remark 1. The time scales T considered in this work have a maximum b and, 
by definition, a(b) = b. For example, let [a, b] — [1, 5] and To = N: in this case 
T = [1, 5] n N = {1, 2, 3, 4, 5} and one has a(t) = t + l,t<= T\{5}, <r(5) = 5. 

Following [3, pp. 2 and 11], we define T fe = T\(p(b),b], T fe2 = (T k ) k and, 
more generally, T fe = ^T fe , for neff. The following standard notation 

is used for a (and p): a°(t) = t, a n (t) = (a o cr™" 1 )^), n G N. 
The graininess function p : T — > [0, oo) is defined by 

= a(t) - t, for all t G T. 

We say that a function / : T — > M is rfeZto differ entiable at t G T fe if there 
is a number f A (t) such that for all e > there exists a neighborhood U of t 
(i.e. £/ = (i - 5, t + 5) (~1 T for some <5 > 0) such that 

\f(a(t)) - f(s) - f A (t)(a(t) - s)\ < s\a(t) - s\, for all s G U. 

We call f A (t) the delta derivative of / at t. 

If / is continuous at t and t is right-scattered, then (see Theorem 1.16 (ii) 
of [3]) 

p(t) 

Now, we define the r th — delta derivative (r G N) of / to be the function 
j/T . T k r _^ Rj provided fA*- 1 ig delta differentiable on T fe \ 

For delta differentiable functions / and g, the next formulas hold: 

r(t) = f(t) + p(t)f A (t), 

(fg) A {t) = f A {t)g°{t)+f{t)g A {t) (3) 

= f A (t)g(t) + r(t)g A (t), 

where we abbreviate here and throughout / o a by f a . We will also write f A 
as f A<y and all the possible combinations of exponents of a and A will be 
clear from the context. 

The following lemma will be useful for our purposes. 

Lemma 1. Let t G H* (t ^ minTJ satisfy the property p(t) = t< cr(t). Then, 
the jump operator a is not delta differentiable at t. 

Proof. We begin to prove that lim s ^ t - a(s) = t. Let e > and take S = e. 
Then for all s G (t — 5, t) we have \a(s) — t\ < \s — t\ < S = e. Since cr(t) > t, 
this implies that a is not continuous at t, hence not delta-diffcrcntiable by 
Theorem 1.16 (i) of [3]. □ 
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A function / : T — ► M. is called rd-continuous if it is continuous in right- 
dense points and if its left-sided limit exists in left-dense points. We denote 
the set of all rd-continuous functions by C r d and the set of all differentiable 
functions with rd-continuous derivative by Cj d . 

It is known that rd-continuous functions possess an antiderivative, i.e. 
there exists a function F with F A = /, and in this case an integral is defined 
by f* f(t)At = F(b) - F(a). It satisfies 

f{r)Ar = ^t)f{t). (4) 

We now present the integration by parts formulas for the delta integral: 

Lemma 2. (Theorem 1.77 (v) and (vi) of [3]) Ifa,b£T and f,g EC\ d , then 

1. S h a f{o{t))g*{t)At = [{fg){t)]t b a j"f A (t)g(t)At; 
2- f a f(t)g A (t)At = [{fg)(t)]\Z b a ~ f a f A {t)g{a{t))At. 

The main result of the calculus of variations on time scales for problem 
(1) is given by the following necessary optimality condition. 

Theorem 1. ([2]) If y* E C\ d is a weak local minimum of the problem 

,6 

C[y{-)} = / L{t,y°{t),y A {t))At ^ min, y(a) = y a , y(b) = y b , 

J a 

then the Euler- Lagrange equation 

L A *{t, yi{t), y A (t)) = y* (*), y*m 

teT k \ holds. 

Remark 2. In Theorem 1, and in what follows, the notation conforms to that 
used in [2]. Expression L A A denotes the A derivative of a composition. 

We will assume from now on that the time scale T has sufficiently many 
points in order for all the calculations to make sense (with respect to this, we 
remark that Theorem 1 makes only sense if we are assuming a time scale T 
with at least three points). Further, we consider time scales such that: 

(H) a{t) = a\t + a for some a\ £ R + and a E M, t 6 [a, b). 

Under hypothesis (H) we have, among others, the differential calculus 
(To — M, ai — 1, ao — 0), the difference calculus (To = Z, a\ = ao = 1) 
and the quantum calculus (To = {q k : k E No}, with q > 1, a\ = q, ao = 0). 

Remark 3. From assumption (H) it follows by Lemma 1 that it is not possible 
to have points which are simultaneously left-dense and right-scattered. Also 
points that are simultaneously left-scattered and right-dense do not occur, 
since a is strictly increasing. 
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Lemma 3. Under hypothesis (H), if f is a two times delta differentiable func- 
tion, then the next formula holds: 

r A (t) = ai f A "(t), t€T k \ (5) 

Proof. We have f aA (t) = [f(t) + p(t)f A (t)] A by formula (3). By the hy- 
pothesis on a, p is delta differentiable, hence [f(t) + p(t)f A (t)] A = f A {t) + 
p A (t)f Aa (t) + n(t)f A2 (t) and applying again formula (3) we obtain f A (t) + 
p A (t)f A °(t)+p(t)f A \t) = f A °{t)+p A {t)f A °{t) = (l + » A (t))f A °(t). Now 
we only need to observe that p A (t) = cr A {t) — 1 and the result follows. □ 

3 Main Results 

Assume that the Lagrangian L(t, u , u\, . . . , u r ) of problem (P) has (stan- 
dard) partial derivatives with respect to uq,...,u t , r > 1, and partial delta 
derivative with respect to t of order r + 1. Let y 6 C 2r , where 

C 2r = {y : T -> R : y A ' r is continuous on T fe2 " } . 

We say that y % e C 2r is a weak local minimum for (P) provided there 
exists 6 > such that C(y*) < C(y) for all y s C 2r satisfying the constraints 
in (P) and \\y - 2/*|| r ,oo < 8, where 

|l?/||r,oo • ^ ^ \\y^ ^||oo) 
i=0 

with j/W = y a%Ar ' and 1 1 2/ 1 1 ^ := sup teT *r- \y(t)\. 

Definition 1. We say that n G C 2r is an admissible variation for problem 
(P) if 

V (a)=0, V (p r -\b))=0 

V Ar ~\a)=0, V Ar - 1 (p r -\b))=0. 
For simplicity of presentation, from now on we fix r = 3. 

Lemma 4. Suppose that f is defined on [a, p 6 (b)] and is continuous. Then, 

under hypothesis (H), jf (b) f(t)n a \t)At = for every admissible variation 
n if and only if f(t) = for all t e [a, p 6 (b)]. 
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Proof. If f(t) = 0, then the result is obvious. 

Now, suppose without loss of generality that there exists to £ [a,p 6 (b)] 
such that /(to) > 0. First we consider the case in which t is right-dense, 
hence left-dense or to = a (see Remark 3). If to = a, then by the continuity of 
/ at t there exists a 5 > such that for all t e [to, to + S) we have f(t) > 0. 
Let us define i] by 

( (t - t ) 8 (t - t - Sf if t E [t , t + (5); 
[0 otherwise. 

Clearly 77 is a C 6 function and satisfy the requirements of an admissible vari- 
ation. But with this definition for 77 we get the contradiction 



L 



pHb) 3 rto+S 

f(t)rf (t)At = / f(t)rf (t)At > 0. 

J to 



Now, consider the case where to 7^ a. Again, the continuity of / ensures the 
existence of a 5 > such that for all t & (to — 5, to + 5) we have f(t) > 0. 
Defining rj by 

= f (t - t + 8) s (t -t - S) s if t e (to - 6, to + 5): 
' otherwise, 



and noting that it satisfy the properties of an admissible variation, we obtain 

fP 5 (6) rt +t 
't - 



/ f(t)^ (t)At = / /(t)rf' (t)At > 0, 

J a Jto—5 



which is again a contradiction. 

Assume now that to is right-scattered. In view of Remark 3, all the points 
t such that t > to must be isolated. So, define 77 such that rf (to) = 1 and 
is zero elsewhere. It is easy to see that 77 satisfies all the requirements of an 
admissible variation. Further, using formula (4) 

/ /(t)rf 3 (t)At = / /(t)Tf' ' (t)At = n(to)f(to)v a3 (t ) > 0, 

J a J to 

which is a contradiction. □ 

Theorem 2. Let the Lagrangian L(t, u , U\, u 2 , u 3 ) satisfy the conditions in 
the beginning of the section. On a time scale T satisfying (H), if y* is a weak 
local minimum for the problem of minimizing 

/ L(t,y cr (t),y a (t),y a (t),y A (t)) At 

J a. 



subject to 
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y(a) = y a , y(p 2 (b)) =yb, 

y A (a) = yl,y A ( P 2 {b))=yl 
y* 2 (a)=yly A2 (p\b)) = yl 
then y % satisfies the Euler- Lagrange equation 

Lu Q {-) ~ <(•) + ^Lf 2 {.) - 1< 3 (.) = 0, t € [a,p e (b)], 
where (•) = (t, yf (t), yf A (t), y°* (t), yf (*)). 

Proof. Suppose that is a weak local minimum of C Let n G C 6 be an 
admissible variation, i.e. n is an arbitrary function such that r], r\ A and r\ A 
vanish at t = a and t = p 2 (b). Define function : E — > R by <P(e) = C{y* +en). 
This function has a minimum at e = 0, so we must have (see [2, Theorem 3.2]) 

<P'(0) = 0. (6) 

Differentiating <P under the integral sign (we can do this in virtue of the 
conditions we imposed on L) with respect to e and setting e — 0, we obtain 
from (6) that 

0=/ {L U0 {-)n° 3 (t)+L Ul (.)^ A {t) 

J a 

+ L U2 {-)n° A \t) + L U3 {-)n A \t))At. (7) 

Since we will delta differentiate L Ui ,i = 1, 2, 3, we rewrite (7) in the following 
form: 

0= / {L U0 W\t)+L Ul (-)v a2A (t) 

J a 

+ L U2 {.)r 1 ° A \t)+L U3 {.)n A \t)}At 
+ n{p\b)){L UQ n° 3 +L UlV ° 2A +L U2 r 1 ° A2 +L U3 n A3 }(p 3 (b)). (8) 
Integrating (8) by parts gives 

0= / {L U0 (.)^ 3 (t)-L A (.)^ 3 (t) 

- L A 2 {.)^ A %t) - L A {-)n A2 %t)}At 

+ [LuA-w^^y + [Lu 2 w A mt p l (h) + K(V 2 W 

+ p(p 3 (b)) [L U0 n° 3 + L ul r? 2A + L U2 n° A2 + L U3 tj a3 } (p 3 (b)). 



t= P s (b) 

t—a 



(9) 
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Now we show how to simplify (9). We start by evaluating if 2 (a): 

if\a)=rf{a)+ l i(a)rf A {a) 

= 77(a) + n(a)rf(a) + p{a)a lV Aa (a) (10) 
= M(a)ai (^(a) + /j(a)^ 2 (a)) 

= 0, 

where the last term of (10) follows from (5). Now, we calculate rf A {a). By 
(5) we have ry crZi (a) = a\r\ Aa (a) and applying (3) we obtain 

a lV A °(a) - ai (rf(a) + ^a)^ (a)) = 0. 

Now we turn to analyze what happens at t = p 3 (b). It is easy to see that if b is 
left-dense, then the last terms of (9) vanish. So suppose that b is left-scattered. 
Since a is delta differentiable, by Lemma 1 we cannot have points which are 
simultaneously left-dense and right-scattered. Hence, p(b), p 2 (b) and p 3 (b) are 
right-scattered points. Now, by hypothesis r] A {p 2 {b)) = 0, hence we have by 
(2) that 

V(p(b))-V(p 2 (b)) 
M (p 2 (6)) 

But r](p 2 (b)) = 0, hence r](p(b)) = 0. Analogously, we have 

v ( P (b))-o^ -0, 

from what follows that rj A (p(b)) = 0. This last equality implies 77(6) = 0. 
Applying previous expressions to the last terms of (9), we obtain: 

r,°\p\b))=r 1 {p{b))=Q, 
r]° 3 (p 3 (b))=r)(b)=0, 

p(p3(b)) 

= 0. 
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In view of our previous calculations, 

r 2 1 t=p 3 (b) A t_„3, M r ~ -|t=p 3 (&) 

LuA-w wi _ + [LuA-^mt: + \L U3 (-)v A (t)} f _ 

J T — O. L -I I — £L 

+ M (P 3 (&)) {i«„^ +L UlV ° 2A +L U2 r 1 ° A ' 2 +L U3 rj A3 } (p 3 (b)) 

is reduced to 1 

L U3 (p 3 (b))v A2 (p 3 (b)) + p(p 3 (b))L U3 (p 3 (b)) V A3 (p 3 (b)). (11) 
Now note that 

v A2 °( P 3 (b)) = v A2 (p 3 (b)) + p(p 3 (b))v A3 (p 3 (b)) 

and by hypothesis r) A ' (p 3 (b)) = r/^ 2 (p 2 (b)) = 0. Therefore, 

p(p 3 (b)) V A3 (p 3 (b)) = -r ] A2 (p 3 (b)), 
from which follows that (11) must be zero. We have just simplified (9) to 

o= / {l U0 (-)^ 3 W-<(-)^ 3 W 

J a 

L A 2 (-)v aA °(t) L A (-)v A2(7 (t)}At. (12) 

In order to apply again the integration by parts formula, we must first make 
some transformations in rf Aa and r) A CT . By (5) we have 

rfA"{t) = -v° 2A (t) (13) 
ai 

and 

V A2 °(t) = ^V° A2 (t)- (W) 

Hence, (12) becomes 
rP 3 (b) 

0= / {L UQ (-W\t)-L A {-W 3 (t) 

J a 

- - 4<(-)^ 2 (t) W (15) 

By the same reasoning as before, (15) is equivalent to 



1 In what follows there is some abuse of notation: L„ 3 (p 3 (6)) denotes L U3 (-)| t 3,^, 
that is, we substitute t in (•) = (t,yf (t),yf A (t),y° A2 (t),y A3 (t)) by p 3 (b). 
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0= / {L U0 (-)v a W-LZW (t) 

J a 

-^ 2 {.)^\t)-^L^{.)r,^\t)\At 



+ l ,{p\b))\L Ua r,° -L^rf -L^rj 



L ~A „a 2 A _ }_ T A a A 1 
_2 L 'uJ\ 



and integrating by parts we obtain 



= 



{L U0 {.)^\t)-L^{.)^\t) 



a i 



-LtW (t) 



t= P 4 (b) 



,2 u 3 



t=p\b) 



(16) 



ai 



Using analogous arguments to those above, we simplify (16) to 

rP 4 (b) 

/ {L ua {-W\t)-L* i {-W\t) 



Calculations as done before lead us to the final expression 

<-P 5 (b) , 



At = 0, 



which is equivalent to 



n 



L „Jr )- <{■) + ±L%(.) - ^Lf 3 {.) \ rTV)Al = 0. (17) 



1 

a 1 



Applying Lemma 4 to (17), we obtain the Euler-Lagrange equation 



U(-) -^(•) + - L ^ 2 2 (-)-^3 3 (-) = 0, te[a,p 6 (6)]. 



□ 
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Following exactly the same steps in the proofs of Lemma 4 and Theorem 2 
for an arbitrary r e N, one easily obtains the Euler-Lagrange equation for 
problem (P). 

Theorem 3. (Necessary optimality condition for problems of the calculus of 
variations with higher-order delta derivatives) On a time scale T satisfying 
hypothesis (H), ify* is a weak local minimum for problem (P), then y* satisfies 
the Euler-Lagrange equation 

E^n 1 ) 2 <(t, y f(t), y f- lA (t),...,y^ r -\t), y f(t))=0, 

(18) 

te[a,p 2r (b)]. 

(i-i)i 

Remark 4- The factor ( ^j-J in (18) comes from the fact that, after each 
time we apply the integration by parts formula, we commute successively a 
with A using (5) (see formulas (13) and (14)), doing this 3 — 

times for each of the parcels within the integral. 
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